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Abstract 

We analytically compute quadratic mass corrections of order 0(a^ s m 2 / s) to 
the absorptive part of the (non-diagonal) correlator of two axial vector cur- 
rents. This allows us to find the correction of order 0(a z s m 2 / M"^) to T{W — > 



hadrons) as well as similar corrections to T(Z — > hadrons). 
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1 Introduction 



Precision measurements of the total and as well as partial Z decay rates have provided 
one of the the most important and, from the theoretical viewpoint, clean determination 
of the strong coupling constant a s with a present value of a s = 0.1202 ± 0.0033 Q. 
Theoretical ingredients were the knowledge of QCD corrections to order oP s in the limit of 
massless quarks plus charm and bottom quarks effects (see, e.g. and references therein). 
These mass corrections which indeed are relevant at the present level of accuracy have 
been calculated up to the order a s s m 2 /s for the vector and a 2 s m 2 q /s for the axial current 
induced decay rate. In this short note the prediction is extended to include aim 2 / s terms 
for the (non-singlet part) of the axial current induced rate. At the same time results 
are obtained for the non-diagonal current correlator with two different masses - a case of 
relevance e.g. for the W decay rate into charmed and bottom quarks. The same formulae 
can also be applied to a subclass of corrections which enter single top production in the 
Drell-Yan like reaction qq tb far above threshold. 

The calculation is based on an approach introduced in Refs. || ||. Knowledge of 
the polarization function to order a 2 , the appropriate anomalous dimensions at order af, 
combined with the renormalization group equation allows one to predict the corresponding 
logarithmic terms of order and hence the constant terms of the imaginary part. The 
first of these ingredients has been available since some time |5], |6|, [7], || while the anomalous 
dimension can been obtained from Ref. || in a straightforward way. 

In this short note only the theoretical framework and the analytical results are pre- 
sented - numerical studies will presented elsewhere. 



2 Renormalization Group Analysis 

In analogy to the vector case, we take as a starting point the generic vector/axial quark 
current correlator V^J V A which is defined by 

UH A ( q ,m u ,m d ,m,ix,a s ) =i [ dxe^(T[j^ A (x)(j^(0)]) 

(1) 

= g, u nty A (Q 2 ) + q,q„nP /A (Q 2 )). 

with Q 2 = —q 2 , m 2 = Y,f m } an d jY^ A = 97^(75)9'- Here q and q' are just two (generically 
different) quarks with masses m u and m d respectively. Note that the vector and axial 
correlators are related through 

n ^(g, m u , m d , m, /i, a s ) = n^(g, m u , -m d , m, /i, a s ) (2) 

The polarization function Hyj A and the spectral density R v ^ A (s) which in turn governs 
the Z and W decays rate obey the following dispersion relation 

r(i)/^,2x _1 f°° J sR v ! A {s,m u ,m d ,m,yL,a s s 



Whereas R y l A physical quantity is invariant under renormalization group transfor- 
mations, the function {T[j^ A {x){j^ A y{Qi) ]) contains some non-integrable singularities 



1 



in the vicinity of the point x = 0. These cannot be removed by standard quark mass and 
coupling constant renormalizations, but must be subtracted independently As a result 
the relevant renormalization group equation assumes the form |2| 



where 



(4) 



(5) 



Here and below the upper and lower signs give the results for vector and axial vector 
correlators respectively. From the identity (Q) we infer that both anomalous dimensions 
7^ and 7*, being not dependent on any masses, also do not depend on the sign. In what 
follows we will denote 

if = lY alld 1m = 1 V J ■ 

The /^-function and the quark mass anomalous dimension 7 m are defined in the usual way 



d/j/ 



a s (3(a s ) = -^2 Pi 



i>0 



2 d 



dys 



m(/i) = m(/i)7 m (a s 



i>0 V 71 J 



(6) 



(7) 



Their expansion coefficients up to order 0(oP s ) are well known [[10], [11], [12|, [U| and read 
(rif is the number of quark flavours) 



/3 = (ll - |n,) /4, fr= (l02-y 
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/) /16, 



C(3) 
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(8) 
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Another useful and closely related object is the correlator of the (pseudo) scalar quark 
currents 



n s / p (Q ; 



,m u ,m d ,m, n,a s 



e l(lx (0\ T [j s/P (x)(j s/p)t ](0) |0). 



(10) 



Scalar and pseudoscalar current correlators are also related in a simple manner: 

n s (Q 2 , m u , m d , m, n) = n p (Q 2 , m u , -m d , m, /x). (11) 

For vanishing quark masses scalar and pseudoscalar correlators are, therefore, identical: 
II = II P and meet the following RG equation 



1 d + 2 7m (a s )V/ p = Q 2 7 f(a s ) ' 



dy? 



16tt 2 



(12) 



The (axial) vector and (pseudo) scalar correlators are connected through a Ward iden- 
tity0 

qtfjtfl* = (m u T m d ) 2 n s/p + (m u T m d )((^ q ) T v0 q 'Vv» > ( 13 ) 

where the vacuum expectation values on the r.h.s. are understood within the framework 
of perturbation theory and the minimal subtractions. Equation (filf ) leads to the following 
relation between the corresponding anomalous dimensions Q|: 

ll V - -7f • (14) 

This relation was used in Ref. Q in order to find the anomalous dimension 7 AA at the a 2 s 
order starting from the results of Ref. || . 

In what follows we will be interested in quadratic mass corrections to the polarization 
operator 11^ which is convenient to represent in the form (m = {m u ,md,m}): 

IT^Q 2 , m, fA, a .) = J^-X/aA^ «•) + ^ U v/a,2(^ m ' «•) + °^)- ( 15 ) 

Here the first term on the rhs corresponds to the massless limit while the second term 
stands for quadratic mass corrections. Note that Tiyj A 2 is a second order polynomial in 
quark masses: a logarithmic dependence on quark masses may appear starting from m 4 
terms onlyQ. 



From the RG equation (Q) we arrive at the following equation for IIyL 2 : 

^ 2 J^v}a,2 = \( m u =F %) 2 7r(« s ) (16) 

2 

or, equivalently, (L q — In 

8 -t(1) _ l /„ ™ \„W (o„. 9 , «„, ^rr(i) 



dT q U v/A,2 = g ( m « T m d )^ v - + 2 lm J U^ /A>2 . (17) 

The last relation explicitly demonstrates that R^ A — the absorptive part of U.yj A 2 - 
depends in order a" on the very function IlyL 2 which is multiplied by at least one factor 
of a s . This means that one needs to know U.yj A 2 up to order a" -1 only to unambiguously 
reconstruct all Q-dependent terms in IIyh 2 to a™, provided, of course, the beta function 
and anomalous dimensions 7 m and 7^ are known to a™. 

This observation was made first in where it was used to find the absorptive part R\ 
in order for the case of the diagonal vector current (that is for the case of m„ = mj. 
In the present paper we will use the results of a recent calculation of ■y^ s || to order a\ to 

determine the absorptive part R^ A to the same order in the general case of non-diagonal 
currents. 



1 Provided of course that one uses a mass independent renormalization scheme like the MS-scheme 
employed in this work. 
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3 Calculation and results 



The result for the function HyJ A 2 in the general non-diagonal case to order a 2 was first 
published in Ref . || . On the other hand, the Ward identity (p~3[ ) expresses the combination 
/Q 2 — IlyL 2 in terms of the massless polarization operator U s known from Refs. |5|, 

(i) 



n 



(i) 

V/A,2 



[?J. A sum of these two functions leads us to the following result for U. V2 
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and m + = m n + m^, Q 2 



Q 2 



-g 2 , all masses as well as QCD 



coupling constant a s are understood to be taken at a generic value of the t' Hooft mass 
//. All correlators are renormalized within MS-scheme. We have also checked ( |i~8|) by a 
direct calculation with the help of the program MINCER [16| written for the symbolic 
manipulation system FORM [|17| . In a particular case of m u = Eq. ([T^) is in agreement 
with Refs. p), g. 

ss 



Now, as was shown in 
the results of IPI we have: 



the anomalous dimension 7^; == — 7? , and, thus, from 
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At last, integrating eq. (O) we find the spectral density R\ in general case to order 
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where the functions r v are 
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4/3 i o J „ oa J 
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The expressions for the hadronic decay rates of the intermediate bosons read: 

T(Z -> hadrons) = rf [ 

+ 



(22) 
(23) 
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T(W -> hadrons) = [ 



2[R (s)+R%(s,0,0,\/m 2 + m 2 c )) 
1 



(25) 



with r r 
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g^p 1 , g I v = Ii - 2Q f sin 2 6 W , g\ = I{ and being the CKM matrix. 
Here Rq(s) is the (non-singlet part) of the ratio R(s) in massless QCD; it was computed 
to in |18|, [L9| and confirmed in |f20|| ; it reads: 
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In deriving Eqs. ( p3j) and (|25|) we have assumed that (i) the top quark is completely 
decoupled (the power suppressed corrections to this approximation start from the order 
-t-o^ and have been studied in Refs. [ 22], p3| , p4|); (ii) all other quarks except for the 



charmed and bottom ones are massless. Note that for the case of diagonal currents there 
exist also so-called singlet contributions to R(s). We will ignore these contributions in 
what follows as they are absent for the case of non-diagonal currents relevant for the 
W-dec&y (a detailed discussion of the Z-decay rate including singlet contributions can be 
found in M). 

Taking into account the peculiar structure of the general result (p0|), the last formula 
can be written in a simpler form, viz. 



T(W -> hadrons) = T^[ 2 (#q(s) + i#(s, 0, 0, ^jm 2 + m 2 )) 

+ R%(s,m ef f, 0,0)]. 



(27) 



Here 



m 



eff= E IVylK + m}) 
i — u,c 
j = d,s,b 

and we have taken into account the fact that 

R v {s, m h rrij, 0) + R A {s, m h rrij, 0) = 2R v {s, Jm 2 + m 2 , 0, 0) = 2R A {s, Jm 2 + m% 0, 0) 
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As a direct consequence of Eqs. (2,|20|) we obtain the following expressions for particular 
functions entering into 
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At last, with rif = 5 and n 2 = s the above formulas are simplified to 
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or, in the numerical form, 
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